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^ ; Abstract 

We estimate the size of the twist-3 effects on deeply virtual Compton scat- 
£f) . tering (DVCS) observables, in the Wandzura-Wilczek approximation. We 

present results in the valence region for the DVCS cross sections, charge asym- 
metries and single spin asymmetries, to twist-3 accuracy. 
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Deeply virtual Compton scattering (DVCS) |}^§ has been studied intensively in recent 
years as a new hard probe to access generalized parton distributions of the nucleon, the 
so-called skewed parton distributions (SPDs). 

After the first experimental evidence for DVCS on the proton [§-0], it is the right time 
to address the problem of the extraction of the twist-2 SPDs from DVCS observables. In 
order to be able to extract the twist-2 SPDs from DVCS observables at accessible values of 
the 'hard' scale Q, one first needs to estimate the effects of higher twist (power suppressed) 
contributions to those observables. The first power correction to the DVCS amplitude is 
of order 0(1/Q), hence it is called twist-3. The twist-3 corrections to the DVCS ampli- 
tude have been derived recently by several groups [§HTl| using different approaches. It was 
shown that the twist-3 part of the amplitude depends on a new type of SPDs which can be 
related to the twist-2 SPDs with the help of Wandzura-Wilczek (WW) relations. These re- 
lations are based on an assumption that the nucleon matrix elements of mixed quark-gluon 
operators are small relative to matrix elements of symmetric quark operators, for details 
H],^ J^O] • Originally, the WW relations were derived for the polarized structure 
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functions g\{x) and gr{x) Recent experimental measurements [15| of these polarized 



structure functions indicate that the WW relations are satisfied to good accuracy. In the 
theory of the instanton vacuum, the smallness of the corrections to the WW relations for 
gr{x) and hi{x) can be related to smallness of the packing fraction of the instantons in the 
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QCD vacuum, for details see fll6|Ji"7 . 

Our aim here is to study the twist-3 effects on DVCS observables for a nucleon target, 
in the Wandzura-Wilczek approximation *. Although on the one hand, the twist-3 effects 
are formally a correction when one aims to extract twist-2 SPDs from DVCS observables, 
they provide on the other hand an additional handle in the extraction of those same twist-2 
SPDs. Indeed, in the WW approximation the twist-3 SPDs are fixed completely in terms of 
twist-2 SPDs and their derivatives. 

In this work, we restrict our analysis to the leading order in a s (Q 2 ), because the NLO 
corrections were computed only for the twist-2 amplitude. Since it was shown that the NLO 
correction to the twist-2 DVCS amplitude can be sizeable |T9|-|2~T|| , the NLO analysis for the 
twist-3 amplitude remains to be investigated. 

In the present work, we firstly give in section the expressions for the DVCS amplitude 
on the nucleon to twist-3 accuracy, as calculated in the WW approximation. 

In section [Til], we describe the modelization of the SPDs that will be used to provide es- 
timates for DVCS observables. We propose a two-component parametrization for the SPDs, 
so as to satisfy the polynomiality conditions. 

In section |V], we show our results for DVCS observables, and discuss the importance of 
the twist-3 effects on the DVCS cross sections, charge asymmetries and single spin asymme- 
tries in the valence region. 

Finally, we give our conclusions in section [V]. 



II. DVCS AMPLITUDE ON THE NUCLEON AT THE TWIST-3 
ACCURACY IN WANDZURA-WILCZEK APPROXIMATION 



The amplitude of the virtual Compton scattering process 

l*(q) + N(p)^ 1 (q') + N(p'), (1) 

is defined in terms of the nucleon matrix element of the T-product of two electromagnetic 
currents : 

T^ = iJ d'x e-^(p'\T [J£ m .(x)r e . m .(0)] \p) , (2) 

where the four-vector index // (u) refers to the virtual (real) photon. 

In the Bjorken limit, where — q 2 = Q 2 — > oo, 2(p • q) — > oo, with xb = Q 2 /2(p • q) 
constant, and t = (p — p') 2 <C Q 2 , the DVCS amplitude on the nucleon to the order 0{1/Q) 
has the form @|TO| : 




*Recently the twist-3 corrections to the DVCS amplitude on the pion target were studied in 
Ref. @ 
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where to the twist-3 accuracy: 
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where p, p' are the momenta of the initial and final nucleon and q, q' are the momenta of 
the initial and final photon respectively. The light-like vectors n and n* are normalized as 
(n ■ n*) = 1. We also introduce the metric and totally antisymmetric tensors in the two 
dimensional transverse plane (£0123 = +1) : 



{-gDjL = -gT + nV + n v n*^ 
The leading order coefficient functions are: 
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In the expression for the DVCS amplitude to the twist-3 accuracy the first two terms 
correspond to the scattering of transversely polarized virtual photons. This part of the 
amplitude depends only on twist- 2 SPDs H, E and H, E and was anticipated in refs. |22| , |23l . 
The third term in Eq. (|3|) corresponds to the contribution of the longitudinal polarization of 
the virtual photon. This term depends only on new 'transverse' SPDs F± and F±. Defining 
the polarization vector of the virtual photon as 



4£ 



(6) 



we can easily calculate the DVCS amplitude for longitudinal polarization of the virtual 
photon (L — > T transition), which is purely of twist-3 : 



L rpfiu 



2 A 

Q J-i 



dx ( Fl C + (x,0-iefF ±k C"(x,0 



(7) 



The skewed parton distributions F^ and F^ can be related to the twist- 2 SPD's H, E, H 
and E with help of Wandzura-Wilczek relations. To derive these relations one assumes that 
the non-forward nucleon matrix elements of gauge invariant operators of the type ipGip are 
small. The WW relations for the case of nucleon SPDs have the form JT0|JT3|] : 
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The functions W±(x,u,£) are the so-called Wandzura-Wilczek kernels introduced in 
Ref. [fL3| . They are defined as : 



w±(x, u, £) = -\e(x> 0- — -r- - e ( x < 0- — 
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The sandwiching between nucleon Dirac spinors is denoted by ((. . .)) = U(p') . . .U(p) and 
M denotes the nucleon mass. The quarks flavor dependence in the amplitude can be easily 
restored by the substitution : 



(13) 



q=u,d,s, 



The amplitude (|3|) is electromagnetically gauge invariant, i.e. 

q^ = {q-A) v T^ = 0, 



(14) 



formally to the accuracy 1/Q 2 - In order to have 'absolute' transversality of the amplitude we 
keep in the expression (0) terms of the A 2 /Q 2 order, applying the prescription of p2| , [2"3f : 
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for the twist-3 terms in the amplitude. Formally such terms are beyond our accuracy and 
they do not form a complete set of 1/Q 2 contributions, but we prefer to work with the DVCS 
amplitude, satisfying Eq. (|15|) exactly. 

The convolution of the leading order Wilson coefficients (||) with the WW kernels can be 
performed analytically (see also [|18j). For example the part of the amplitude which contains 
WW kernels can be simplified using the following result for the integration over x : 
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Note that all integrals over u are convergent and the strongest (integrable) singularity of 
the integrand is logarithmic only. 



III. MODELIZATION OF SPDS 

In order to interpret DVCS observables, we next discuss the parametrization of the SPDs. 
As we address here hard exclusive reactions in the valence region at relatively small values 
of —t, we factorize the t-dependence of the SPDs by the corresponding form factor, so as 
to satisfy the first sum rule. We stress that the factorization of the t-dependence of SPDs 
is not expected to hold when increasing the value of t. For example, the calculations of the 
SPDs in the chiral quark soliton model PH| showed that at t about —0.8 GeV 2 the deviations 



from the factorization ansatz for the t-dependence of SPDs is as large as about 25% in some 
x-range. In the calculations presented here, we limit ourselves to relatively small values 
of —t = 0.25 GeV 2 and assume the factorization ansatz. A more realistic parametrization 
of SPDs at larger values of t, remains to be done in a future work. We next discuss the 
resulting t-independent part. 

For the function H q (for each flavor q), the t- independent part is parametrized by a 
two- component form : 

H q (x,0 = H* DD (x,t) + ^rD(^), (17) 
where H q DD is the part of the SPD which is obtained as a one-dimensional section of a 



two-variable double distribution F q [24], imposing a particular skewedness £, as 



H q DD (x,£)= [ d(3 f m da5(x-(3-aO Fi(fra) , (18) 
J-i •/ -1+101 
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and where the D-term contribution^] D completes the parametrization of SPDs, restoring 



the correct polynomiality properties |25[ of SPDs ||26|| . The D-term contribution to SPDs 
has a support only for \x\ < |£|, so that it is 'invisible' in the forward limit. The D-term is 
an isoscalar contribution, and adds therefore the same function for each flavor (in Eq. (|17|), 
Nf = 3 is the number of active flavors). 

For the double distributions, entering Eq. (|18D, we use the following model suggested in 
El : 



(A a) ~ 2^(6+1) (l-|^|)2M-i M) > ( 19 ) 

with q{(3) is the corresponding forward parton distribution ensuring the correct forward limit 
for the SPD H(x,0) = q(x). We use the phenomenological forward parton distributions 
(including their evolution) as input. In Eq. (|19|), the parameter b characterizes the strength 
of the £ dependence of the SPD H(x,£), the limiting case b — > oo corresponds to the £ 
independent SPD H(x,£) = q(x). The power b in eq. ( |I9"D is a free parameter for the valence 
contribution (b va i) and for the sea/antiquark contribution (b sea ) to the SPD, which can be 
used as fit parameter in this approach to extract SPDs from DVCS observables. The twist-2 



DVCS predictions of Ref. [pq| , correspond to the choice : b va i = b sea = 1.0. In the following, 
we also show all predictions for the value : b va [ = b sea = 1.0, in order to have a point of 
comparison, and refer to a future work for a more systematic study of the dependence of 
DVCS observables on the shape of the profile function. 

The D-term contribution to the singlet SPDs H(x,£) and E(x,£) has the following form : 

rD-term' ^ ~ ' X 



H D -^(x,0 = e^-\x\)D[-) , (20) 



S D " to ( a ; ) e) = -e(e-N)Df|l . (21) 

Note that the D-term contribution is cancelled in the combination H + E to ensure the 
polynomiality condition (for odd N) : 



J dx x N (H(x,£) + E(x, £)) = polynomial of the order (N — 1) in £ . 



(22) 



Note that the iV-th Mellin moments of H and E separately are polynomials of order N + 1 
in £, only in the combination ( |2"2"D the highest powers are cancelled. As model for the SPD 
E(x,£), we shall only use the contribution of the D-term (|21|). More detailed studies of 
DVCS observables with the function E(x,£) computed in the chiral quark-soliton model 
29| , |30[] will be published elsewhere. 



The function D(a) can be expanded in odd Gegenbauer polynomials : 



TThe name "D-term" is derived from (more or less arbitrary) notations used in Ref. [26 and 
to some extent is confusing because of its similarity to the terminology used in supersymmetric 
theories, especially when the first term in Eq. (117]) is called "F-term". 
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(23) 



For the moments di,d^,d 5 , we use the estimate which is based on the calculation of SPDs in 
the chiral quark soliton model [E9| at a low normalization point fi « 0.6 GeV, which gives : 



d x « -4.0, 



d, « -1.2, 



4 w -0.4, 



(24) 



and higher moments (denoted by the ellipses in Eq. fl23|) ) are small and are neglected in 
the following. Notice the negative sign of the Gegenbauer coefficients for the D-term in 
Eq. (^4[), as obtained in the chiral quark-soliton model. 

It is interesting to note that for the case of SPDs in the pion, the value of the coefficient 
d\ in the parametrization of the D-term (|23| ) can be computed in a model independent way. 
To do this, we can use the soft-pion theorem for the singlet SPD in the pion derived in 
This soft-pion theorem states that the singlet SPD in the pion vanishes for £ = ±1 : 



^#W(x,£ = ±l,t = 0) 







(25) 



Hence (see also p6|), 



J X dxx ^#W(z,£,f = 0)j = (l-£ 2 ) Mi 



(26) 



Evaluating Eq. ([26]) at £ = 0, determines : 



i 

M 2 V = I dxx V[g W (a;) + q {n) (x)] 



(27) 



being the fraction of the momentum carried by the quarks and antiquarks in the pion. As 
the highest power in £ in Eq. (|26"D (i.e. the term in £ 2 ) originates solely from the D-term, 
one easily obtains the following expression for pion D-term : 



5M 2 Q 
4 



(1-**) Cf /2 (z) + 



(28) 



We see therefore that the first Gegenbauer coefficient of the pion D-term is negative and 
strictly nonzero. The coincidence of the signs in the nucleon and pion D-terms hints that 
the D-term in the nucleon is intimately related to the spontaneous breaking of the chiral 
symmetry. 

When using the parametrization of Eq. fll?D for the pion, the soft pion theorem (|25"D 
imposes a condition between both terms, and fixes the pion D-term completely in terms of 
the isosinglet double distribution as follows : 



4n what follows we neglect the scale dependence of the result (^4|) as the uncertainties in modeling 
of the D-term are larger than the logarithmic scale dependence, for the scales considered in our 
calculations. 
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pi fl-\p\ 

D ( ^{z) = - dp / da 6(z-p-a) VF 9 (/3,a) . (29) 

J-l J-l+W g 

Therefore, when studying DVCS off the pion it is enough to model the double distribution, 
the corresponding D-term is restored by using the soft pion theorems. For the nucleon, both 
terms in Eq. ( |TTD are parametrized independently. 

For the calculation of the twist-3 terms, we need derivatives of SPDs H(x,£) which can 
be computed in terms of derivatives of double distributions, using the obvious result : 



{x^ + ^)H DD (x,0 = ~H DD (x,0 + f*Pj da 8(x-P- afl dF( ^ a) . (30) 



For the model of Eq. (Jl9|), the derivative with respect to a in Eq. fl3"0|) is easily performed: 
dF(p,a) ba T(26 + 2) [(1 - \(3\) 2 - a 2 }'' 1 



da 2 2b r 2 (6 + l) (1- \P\) 2b+l 



q(P) . (31) 



Note that although generically the derivatives of a SPD H(x,£) with respect to x and £ are 
discontinuous at the points x = ±£, the combination (xd x + £d^)H DD (x, £) is continuous in 
these points. The contribution of the D-term to the 'transverse' SPDs F± and F± has been 



computed in Ref. [13] with the result : 



F ww D- tcrm(a , = -J^((i))d (j) e(\x\ < o , 

F M WW D " term (x,0 = 0. (32) 

Therefore the contribution of the D-term to the DVCS amplitude can be easily added sep- 
arately from the contribution of the double distribution part Hdd of Eq. ([0J). Obviously 
the D-term contributes to the real part of the amplitude only. 

In Fig. H| we plot for illustrative purposes H u (x,£) obtained from the model (|T9"D with 
frvai = &sca = 1) using the MRST98 parametrization |yj for the forward quark distribution 



as input. We also show separately the contribution of the D-term (23) to this function, 
as well as the combination (1 + xd x + C,d^)H u (x,^) which enters the twist-3 SPDs in WW 
approximation. 

To model the SPD H, we make use of its representation in terms of double distributions 
in the same way as for the function H, see Eqs. (0,|18l[T9D. For the forward polarized parton 



distributions, we use (as in ||28|| ) the parametrization of Ref. [32]. Note that calculations in 
the chiral quark-soliton model of the functions H |36[] show that the corresponding D-term 
is very small and can be neglected. 



As model for E we take the contribution of the pion pole |p3|-36j of the form 



E*» polc (x,0 = ( f) 0{\x\ < , (33) 

where g& is the axial charge of the nucleon and (p n (z) is the pion distribution amplitude 
(DA). The calculations are performed with an asymptotic DA for the pion. Again the 
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contribution of the pion pole to the 'transverse' functions F± and F±_ has been computed in 
13fl with the result : 



F W,pion P°le( X) £) = 0> 

pWW^ pole (X;0 = ^ H)) ^L 1^ ^ < . (34) 



IV. RESULTS AND DISCUSSION 

In this section, we present our results for the DVCS observables. We give all results for 
the invariant cross section of the ep — > epj reaction, which is differential with respect to Q 2 , 
Xb, t, and out-of-plane angle $ ($ = 0° corresponds to the situation where the real photon 
is emitted in the same half plane as the leptons). The invariant ep — > epy cross section is 
given by : 



<Ia I x u y 2 f i AM 2 x\ ^ " ' 2 



dQ 2 dx B dtd$ (2tt) 4 32 Q 4 \ Q 



Tbh + T FV cs 



(35) 



where M is the nucleon mass, y = (p ■ q)/{p- k), and k is the initial lepton four-momentum. 
In the ep —>■ epj reaction, the final photon can be emitted either by the proton or by the 
lepton. The former process is referred to as the fully VCS process (amplitude T FV cs in 
Eq. ([35])), which includes the leptonic current. The process where the photon is emitted 
from the initial or final lepton is referred to as the Bethe-Heitler (BH) process (amplitude 
Tbh in Eq. (]35f)), and can be calculated exactly. For further technical details of the ep — > epj 
reaction and observables, we refer to Ref. ^2] . 



When calculating twist-3 effects in DVCS observables, we show all results with the exact 
expression for the BH amplitude, i.e. we do not expand the BH amplitude in powers 1/Q. 
Also in the actual calculations, we use the exact kinematics for the four-momenta of the 
participating particles (the kinematics in Eq. (H) were shown to twist-3 accuracy only for 
simplicity of the presentation). In this way we take (partially) kinematical higher twists 
into account. Furthermore, when referring to the twist- 2 DVCS results, we include those 
higher twist effects which restore exact transversality of the amplitude as expressed through 
Eq. (pl|). Similarly, when referring to the twist-3 DVCS results (providing from a longitudi- 
nally polarized virtual photon), we include the gauge restoring higher twist terms. 

In Fig. [| we show the ^-dependence of the ep — > epry cross section for a lepton (either 
electron or positron) of 27 GeV (accessible at HERMES). The kinematics corresponds to 
the valence region (xb = 0.3) and to a ratio t/Q 2 = 0.1 (Remark that increasing the ratio 
t/Q 2 , increases the higher twist effects). 

It is firstly seen from Fig. |2], that in these kinematics, the pure twist- 2 DVCS process 
without D-term contribution (dashed curves) dominates the ep — ► epj cross section com- 
pared to the BH process (which is sizeable only around $ = 0°, where it reaches its maximal 
value). In absence of the BH, the twist-2 DVCS cross section would give a "^-independent 
cross section which practically saturates the ep — > epr) cross section at $ = 180°, where the 
BH is vanishingly small. The $- dependence of the dashed curve in Fig. |^ can be understood 
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as the sum of this constant twist-2 DVCS cross section, the cross section for the BH process 
(whose amplitude is purely real), and the interference of the BH with the relatively small 
real part of the DVCS amplitude (in the valence kinematics, xb — 0.3, shown in Fig. |2], the 
ratio of real to imaginary part of the DVCS amplitude without D-term, is around 15 %). 

When adding to the twist-2 DVCS amplitude the purely real D-term contribution, the 
resulting cross section for the full twist-2 DVCS process is given by the dashed-dotted curves 
in Fig. 0. At $ = 180° (where the BH 'contamination' is very small), the predominantly 
imaginary DVCS amplitude (in absence of the D-term contribution), and the purely real 
D-term amplitude have only a small interference. In this region, the DVCS cross section is 
enhanced by about 10 %, when using the chiral quark soliton model estimate of Eqs. ( p3| , |24D 
for the D-term. Going to $ = 0°, the purely real D-term amplitude interferes maximally 
with the BH. This interference is destructive for the electron reaction and constructive for 
the positron reaction. Therefore, the effect of the D-term can be very clearly seen in the 
^-dependence of the charge asymmetry as shown on the lower panel of Fig. ^. Including the 
D-term contribution, the full twist-2 DVCS charge asymmetry changes sign and obtains a 
rather large value (~ 0.15) at $ = 0°. On the other hand, at $ = 180°, where the interference 
with the BH (and hence the difference between e~ and e + ) is small, the charge asymmetry 
is correspondingly small. The pronounced ^-dependence of the charge asymmetry and its 
value at $ = 0°, provides therefore a nice observable to study the D-term contribution to 
the SPDs, and to check the chiral quark soliton model estimate of Eqs. (^,^4|) . 

Adding next the twist-3 effects in the WW approximation (full curves), it is seen from 
Fig. H] that they induce an additional (approximate) cos $ structure in the cross section. The 
twist-3 effects would induce an exact cos $ structure only when the BH amplitude is approx- 
imated by its leading term in an expansion in 1/Q, neglecting its additional ^-dependence. 
In all calculations, we keep however the full ^-dependence of the BH, due to the lepton 
propagators, which complicates the interference at the lower Q. One sees from Fig. [| that 
the interference of the twist-3 amplitude with the BH + twist-2 DVCS amplitude is destruc- 
tive for $ < 90°, and constructive for $ > 90°. At $ = 0°, the twist-3 effects reduce the full 
twist-2 DVCS cross section (including the D-term) by about 25 %, whereas at $ = 180°, 
they largely enhance the twist-2 DVCS cross section (by about 55 %). 

To further study the twist-3 effects on the ep epy cross section, we show in Fig. ^ the 
cross section and charge asymmetry at the same E e , xb, and t, as in Fig. 0, but at a value 
Q 2 = 5 GeV 2 . One firstly sees that going to higher Q 2 , at the same E e , xb, and t, enhances 
the relative contribution of the BH process compared to the DVCS process. Consequently, 
the cross section follows much more the $-behavior of the BH process. For the twist-2 cross 
section, one sees again clearly the effect of the D-term which leads, through its interference 
with the BH amplitude, to a charge asymmetry of opposite sign as compared to the one for 
the twist-2 DVCS process without D-term. 

The twist-3 effects in the kinematics of Fig. ^ where to ratio t/Q 2 is only half the value 
of Fig. are correspondingly smaller. It is furthermore seen that the twist-3 effects in- 
duce an (approximate) structure ~ A cos(2$) in the charge asymmetry (A w - 0.08 in the 
kinematics of Fig. |3|). Only for out-of-plane angles $ > 130° one sees a deviation from the 
simple cos(2<3>) structure induced by the twist-3 amplitude, due to the more complicated 
^-dependence when calculating the BH amplitude exactly. 

In Figs. f| and |5], we compare the ^-dependence of the DVCS cross section and of the 
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single spin asymmetry (SSA) corresponding to a polarized electron beam, for the same val- 
ues of Q 2 and t as in Figs. || and [| but for a value of xb = 0.15. Note that, due to parity 
invariance, the SSA is odd in $. We therefore display only half of the $ range, i.e. $ 
between 0° and 180°. 

When comparing Figs. |2|, |3] and Figs. |], |5], it is firstly seen that the cross sections for the 
ep — > epj process increases strongly when decreasing at fixed Q 2 and fixed t, mainly due 
to the growth of the BH amplitude. Due to this large BH amplitude, the relative twist-3 
effects in Figs.^, [5] are smaller than the ones in Figs. |[ However, the large BH amplitude 
leads to a large value for the SSA through its interference with the DVCS process. At (pure) 
twist-2 level, the SSA originates from the interference of the imaginary part of the DVCS 
amplitude and the (real) BH amplitude. In case the BH is approximated by its leading term 
in an expansion in 1/Q, the twist-2 SSA displays a pure sin$ structure. Due to the more 
complicated ^-dependence of the BH at the lower values of Q 2 , this form gets distorted and 
its maximum displaced, as is seen by a comparison of the twist-2 SSA in Fig. ^ (Q 2 = 2.5 
GeV 2 ) and Fig. [5| (Q 2 = 5 GeV 2 ). Note that, for practical considerations, our "twist-2" 
DVCS calculations include kinematical higher twist terms as well as the gauge restoring 
terms according to Eq. (0). Their effect can be seen in the slight change in the SSA (of the 
percent level), due to the DVCS process by itself (i.e. when increasing the real part of the 
DVCS amplitude by adding the D-term contribution, the curves for the SSAs are slightly 
displaced). A more systematic treatment of target mass corrections of order M 2 /Q 2 , and 
corrections of order t/Q 2 for DVCS, still remains to be done. 

We next discuss the twist-3 effects, calculated in WW approximation, on the DVCS SSA. 
One sees from Fig. |4] that the twist-3 corrections induce an (approximate) sin(2$) structure 
in the SSA. The amplitude of the sin(2$) term is however rather small and the twist-3 effects 
change the SSA by less than 5 % in the kinematics corresponding to t/Q 2 = 0.1. It was 
checked that at xb = 0.3 and for a value t/Q 2 = 0.1, the twist-3 effects on the SSA are of 
similar size. One therefore observes that although the twist-3 effects in WW approximation 
can provide a sizeable contribution to the real part of the amplitude (Fig. |2|), they modify 
the imaginary part, and hence the SSA, to a much lesser extent. 

Finally, we show in Fig. [(], the corresponding cross sections at 11 GeV (JLab), for the 
same Q 2 , xb and t as in Fig. [2| From the cross sections, one firstly observes again the rather 
large effect of the D-term, which changes the twist-2 cross section by around 25 %. The 
relative twist-3 effects on the cross section are smaller in this case than in Fig. 0, as they are 
on top of a larger BH amplitude. For the SSA, it is again observed that the twist-3 effects 
induce an approximate sin(2$) structure, with amplitude of less than 5 % in the valence 
kinematics considered here. 

V. CONCLUSIONS 

We have estimated the size of twist-3 effects on DVCS observables in the Wandzura- 
Wilczek (WW) approximation, which allows to express the new twist-3 SPDs in terms of 
twist-2 SPDs and their derivatives. The twist-3 effects in WW approximation display there- 
fore a new sensitivity to the shape of the twist-2 SPDs and provide an additional handle in 
the extraction of twist-2 SPDs from DVCS observables. The WW approximation relies on 
the assumption that nucleon matrix elements of quark-gluon operators are small and can be 
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neglected compared to the matrix elements of quark operators. Recent data on polarized 
structure functions indicate that this approximation holds to good accuracy. Therefore, the 
WW approximation was considered in this paper as a benchmark to study the size of the 
twist-3 effects on DVCS observables. It remains for future study to quantify the remaining 
twist-3 effects due to matrix elements of mixed quark-gluon operators. 

To provide estimates and to interpret DVCS observables, we proposed a two-component 
parametrization for the SPDs of the nucleon as a sum of a double distribution part and a 
D-term, so as to satisfy the polynomiality conditions. The D-term contribution is intimately 
related to the spontaneous breaking of the chiral symmetry, and was estimated using the 
chiral quark soliton model. It was shown that the effect of the D-term, which gives a purely 
real contribution to the DVCS amplitude, can be very clearly seen in the dependence of the 
ep — > ep7 cross section on the out-of-plane angle $. Through its interference with the BH 
process, the D-term contribution was seen to change the sign of the DVCS charge asym- 
metry, and to lead to a large value for the charge asymmetry at $ = 0°. The pronounced 
charge asymmetry provides therefore a clear signature to study the D-term and to check the 
chiral estimate given here. The importance of a quantitative understanding of the D-term is 
especially emphasized by the fact that it contributes considerably to the DVCS cross section 
but drops out in the combination H + E in Ji's angular momentum sum rule 0. 

The twist-3 effects were seen to give a sizeable change of the real part of the DVCS 
amplitude in the valence region, and they introduce an (approximate) cos $ structure in the 
cross section. In kinematics where the DVCS amplitude is comparable to the BH amplitude, 
one gets sizeable effects. In particular, it was shown for the e~p — > e~pry reaction at E e = 
27 GeV, in valence kinematics (xb — 0.3), and for a value t/Q 2 = 0.1 (the value t/Q 2 governs 
kinematically the size of higher twist effects), that the twist-3 effects reduce the full cross 
section around $ = 0° by around 25 %, and largely enhance the cross section at $ = 180° 
(where the BH is very small). On the charge asymmetry, the twist-3 effects introduce an 
(approximate) cos(2$) structure, with amplitude around 10 % in those same kinematics. 

For the single spin asymmetry, which is proportional to the imaginary part of the DVCS 
amplitude, the twist-3 effects introduce an (approximate) sin(2$) structure whose amplitude 
is however rather small, so that the SSA gets changed by less than 5 % in typical valence 
kinematics for a value t/Q 2 = 0.1. 

In summary, we found that the ep — > epry cross section and charge asymmetry in the 
valence region, display a sensitivity on twist-3 effects, in particular through the real part of 
the amplitude. In the WW approximation, this sensitivity to the shape of the twist- 2 SPDs 
may be exploited in fact in the extraction of those same twist-2 SPDs from DVCS observ- 
ables. Such a systematic study of the extraction of twist-2 SPDs from DVCS observables 
to twist-3 accuracy in the WW approximation, seems to be a promising subject for future 
work. 
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FIGURES 
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FIG. 1. SPD H u (x,£,) for the u-quark obtained from the double distribution model ( fl^ j (solid 
curve); (1 + xd x + ^d^)H u (x,^) (dotted curve), and D-term ( pffi ) (dashed- dotted curve) at a fixed 
value of £ = 0.3. 
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FIG. 2. Invariant cross section for the e~p — > e~p7 reaction (upper panel) and DVCS charge 
asymmetry (lower panel) at E e = 27 GeV, for the DVCS kinematics as indicated in the figure. Dot- 
ted curves : BH contribution; dashed curves : BH + twist-2 DVCS (without D-term); dashed-dotted 
curves : BH + twist-2 DVCS (with D-term); full curves : BH + twist-3 DVCS. 
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e7e + + p — > e7e + + p + y 
E p = 27 GeV, Q 2 = 5 GeV 2 , x R = 0.3, t = -0.25 GeV 2 



> 

o 

a 



0.01 



0.008 



O 0.006 



K « 0.004 

CM 

-5 0.002 



T3 




20 



40 



60 80 100 120 140 160 180 



+ 

+ 



to 



0.25 
0.2 

0.15 
0.1 

0.05 


-0.05 
-0.1 

-0.15 
-0.2 

-0.25 




20 40 60 80 100 120 140 160 180 

<t> (deg) 



FIG. 3. Same as Fig. |, but for Q 2 = 5 GeV 2 . 
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FIG. 4. Invariant cross section for the e~p — > e~pj reaction (upper panel) and DVCS SSA 
(lower panel) at E e = 27 GeV, for the DVCS kinematics as indicated in the figure. Curve conven- 
tions as in Fig. H. 
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FIG. 5. Same as Fig. % but for Q 2 = 5 GeV 2 . 
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FIG. 6. Invariant cross section for the e~p — > e~pj reaction (upper panel) and DVCS SSA 
(lower panel) at E e = 11 GeV, for the DVCS kinematics as indicated in the figure. Curve conven- 
tions as in Fig. H. 
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